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Abstract 

The potential of B-physics in ultrarelativistic heavy ion collisions is discussed. The different 
mechanisms of heavy quark production at high energies are considered. We analyze the sensitivity 
of high-mass pairs from BB semileptonic decays and secondary J/^'s from single B decays 

to the medium-induced bottom quark energy loss at LHC energies. 



I. INTRODUCTION 

The production of heavy b- and c-quarks in hadronic interactions at very high energies allows one to 
study the dynamics of hard processes describing the standard as well as new QCD-physics dJ^]. This 
production can be described completely in the frame of perturbative QCD under condition the mass of 
quark Mq is much larger than characteristic QCD confinement scale, Mq ^> Aqcd- The specific interest 
for heavy quark production in ultrarelativistic nuclear interactions is due to the intriguing possibility 
to study behaviour of massive colour charge in super-dense QCD-matter - quark-gluon plasma (QGP), 
the search and investigation of properties of which being one of the goals of modern high energy physics 
(see, for example, reviews |3p|6|]). Experimental data obtained at CERN-SPS and RHIC-BNL can be 
interpreted as a result of QGP formation in heavy ion collisions, although alternative explanations can 
not be fully dismissed 0. It is expected that the most central heavy ion collisions at LHC collider 
might produce practically "ideal" quark-gluon plasma at extremely high energy density up to £o ~ 0.5 
TeV/fm 3 3> e C rit ~ 1 GeV/fm 3 The inclusive cross section for 6-quark production at LHC (y/s — 5.5 A 
TeV for Pb— Pb) will be large enough for systematical studies of different aspects of -B-physics, while 
at RHIC (y/s — 200A GeV for Au— Au) the copious production only of c-quarks can be expected. 
The heavy quark pairs, created at the very beginning of the collision process, propagate through the 
dense matter and interact strongly with constituents of the medium. In-medium gluon radiation and 
collisional energy loss (see review || and references therein) of heavy quarks can result in experimentally 
observed modification of high mass dilepton spectra p0|-[l3| . It was also predicted recently that the finite 
quark mass effects can lead to a relative suppression of medium-induced radiation of heavy quarks with 
corresponding enhancement of heavy-to-light D(B)/n ratio Since the parton rescattering intensity 
strongly increases with temperature, formation of a "hot" QGP at initial temperatures up to To ~ 1 
GeV at LHC || should result in much larger parton energy loss compared to " cold" nuclear matter or a 
hadronic gas. 

II. MECHANISMS OF HEAVY QUARK PRODUCTION AND DIMUON SPECTRA 

The diagrams for different mechanisms of heavy quark production in high energy hadronic interactions 
are presented in figure 1. Following the classification and terminology of paper 0, we can distinguish 
three classes of such processes with different number of heavy quarks in hard sub-processes vertices (2, 
1 or 0): "pair creation" (leading order diagrams are shown on the top of fig.l); "flavour excitation" 
(one heavy quark is produced in vertex of hard process and another quark is created from initial state 
parton shower); "gluon splitting" (both heavy quarks are produced from final state parton shower). At 
present accelerator energies (including RHIC) the bulk of heavy quarks are produced due to direct hard 
scattering ("pair creation"). However, the prediction for LHC is that the contribution of gluon splittings 
in initial- or final-state shower evolution to heavy flavour yield can be a significant, about 90% for the 
whole kinematical range [0. 



*Talk given at International Workshop on Physics of the Quark-Gluon Plasma, Palaiseau, France, September 
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Since heavy quark pairs are produced at the very beginning of the nuclear collisions, its propagate 
through the dense medium. They finally form B and D mesons by "capturing" u, d or s quarks during 
the hadronization stage. These mesons will decay with the average meson lifetimes ct b ± = 496 pm, 
ct b o — 464 nm, ct d ± = 315 pm and ct d o — 124 /im. We note that w 20% of B mesons and ~ 12% of D 
mesons decay to muons. About half of the muons from B decays are produced through an intermediate 
D and contribute to the softer part of the pT-spectrum. 




FIG. 1. Examples of diagrams for different heavy FIG. 2. Initial invariant mass distribution of p + p~ 

quark production mechanisms from top to bottom cor- pairs from BB decays for different production mech- 

respond to the "pair creation", "flavour excitation" and anisms: "pair creation" (solid), "flavour excitation" 

"gluon splitting". (dashed) and "gluon splitting" (dotted). 



In order to estimate dimuon spectra at LHC we used the initial momenta spectra, heavy quark cross 

sections in NN collisions and meson fragmentation scheme from PYTHIA5.7 |1| with the default 
CTEQ2L parton distribution function (PDF) and including initial and final state radiation in vacuum 
which effectively simulates higher-order contributions to heavy quark production. The initial distribution 
of QQ pairs over impact parameter b of A— A collisions (without any nuclear collective effects) is obtained 

by multiplying er^^ by the number of binary nucleon-nucleon sub-collisions: 



'-^.(b,VS) = T AA (b)^ Q N (V~s) 



±T AA (b)*% N (VS) 



(1) 



with the total inelastic non-diffractive nucleon-nucleon cross section is cr^ N — 60 mb at ^/s = 5.5 
TeV. The standard Wood-Saxon nuclear overlap function is T AA (b) — J d 2 sT A (s)T A (\b — s|) where 
T A (s) = A J dzp A {s, z) is the nuclear thickness function with nucleon density distributions p A (s, z). 

Figure 2 shows the initial p + p" invariant high-mass spectra from dominate source - semileprtonic BB 
decays - in the CMS experiment (T^| kinematical acceptance, p^, > 5 GeV/c and |ry M | < 2.4. The total 
impact parameter integrated rates are normalized to the expected number of Pb— Pb events during a two 
week LHC run, R = 1.2 X 10 6 s, assuming luminosity L — 10 27 cm _2 s _1 |l(| to that 

N( f i+ t i-) = Ra% + /~L. 

. We can see the approximately equal contribution from "pair creation" and "flavour excitation" at 
high masses, while the relative contribution from "gluon splitting" being decreased with M increasing. 
The integrated from M p + AI - > 20 GeV/c 2 cross sections are 26.5 pib (3.2 x 10 4 events per two week), 
20.5 pb (2.5 x 10 4 events) and 9.5 pb (1.15 x 10 4 events) for "pair creation", "flavour excitation" and 
" gluon splitting" respectively. In the remainder of the discussion we will consider as signal only dimuons 
from bb decays, because the contribution of cc fragmentation into kinematical region of interest is only 
about 1/5 part of total rate. Moreover, the medium-induced charm quark energy loss can be significantly 
larger than the 6-quark loss due to the mass difference |l4| , resulting in an additional suppression of the 
DD -> p+p- yield. 
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Note that main correlated background here - Drell-Yan prompt dimuons - is unaffected by medium- 
induced final state interactions. These dimuons are directly from the primary nuclear interaction vertex 
while the dimuons from B and D meson decays appear at secondary vertices on some distance from the 
primary vertex. The path length between the primary vertex and secondary vertices are determined by 
the lifetime and 7-factor of the mesons. This fact allows one cut to suppress the Drell-Yan rate by up to 
two orders of magnitude using the dimuon reconstruction algorithm based on the tracker information on 
secondary vertex position [T^ ] . 

The main uncorrelated background - random hadronic decays and muon pairs of mixed origin - appears 
also in the like-sign dimuon mass spectra and can be subtracted from the total fj, + [i~ distribution using 
the and yT fi~ event samples. 

Another process of particular interest is secondary J/ip production [ Ell . The branching ratio B — > 
J/tpX is 1.15%. The J/ip's subsequently decay to dimuons with a 5.9% branching ratio. The main 
correlated background here - primary J/ip production - as in the case of Drell-Yan dimuons, can be 
rejected using vertex detector. The estimated in the same kinematical region cross section for secondary 
J/ip(—> production from "pair creation" is only 17% (ss 10.5 /ib, i.e. 1.3 x 10 4 events per 

two weeks) of total number of secondary J/ip's, while the contribution from "flavour excitation" and 
"gluon splitting" being w 25.5 fjh each. Since we consider here the region of much lower dimuon masses, 
Mp+p- = Mj/^ = 3.1 GeV/c 2 , the contribution of " showering" bb pairs is significantly larger as compared 
with the dimuons from semileptonic BB decays. On the other hand, the spectra shape of secondary J/ip's 
is very similar for different bb production mechanisms, because it carries information about spectra of 
single 6-quarks, which are practically the same for the different sources 0]. The significant difference 
appears only for the quark-antiquark correlations changing variables sensitive to the QQ kinematics 
(diquark invariant mass, azimuthal correlations). 

Note that there are theoretical uncertainties in the bottom and charm production cross sections in 
N— N collisions at LHC: the absolute dimuon rates depend on the choice of PDF, the heavy quark mass, 
the i3-meson fragmentation scheme, next-to-leading order corrections, etc. It is therefore desirable that 
dimuon measurements in pp or dd collisions are made at the same or similar energy per nucleon as in the 
heavy ion runs. 



III. MEDIUM-INDUCED ENERGY LOSS OF BOTTOM QUARKS AND DIMUON SPECTRA 

The details of our model for the heavy quark production and passage through a gluon-dominated 
plasma, created in initial nuclear overlap zone in minimum bias Pb— Pb collisions at LHC, can be found 
in our work |jl3f . We treat the medium as a longitudinally expanding fluid with parton production on a 
hyper-surface of equal proper times t ]17t and perform a Monte-Carlo simulation of the mean free path 
A of heavy quark in scattering-by-scattering scheme. Then the basic kinetic integral equation for total 
energy loss AE as a function of initial energy E and path length L has the form 

L 

r,s f , dP(x),, s dE(x,E) dP(x) 1 , /w „ 

AE(L, E) = / dx— ^A(x) V , — 7^ = tt^ cxp {-x X(x)), (2) 
J dx dx dx X{x) 



where x(t) is the current transverse coordinate of a quark. 

Note that although some attempts have been made to calculate medium-induced heavy quark energy 
loss for quarks of mass M q (see p^ , p^ , |l^ | for discussion) , a full description of the coherent gluon radiation 
from a massive colour charge still lacking. There are two extreme limits for energy loss by gluon radiation. 
In the low px limit, px < M q , medium-induced radiation should be suppressed by the mass, while the 
ultrarelativistic limit, px — > 00, corresponds to the radiation spectrum of massless quarks. In our case, the 
main contribution to high-mass dimuon and secondary charmonium production is due to 6-quarks with 
"intermediate" values of pr > 5 GeV/c, expected to be rather close to the incoherent regime. In order to 
estimate the sensitivity of the dimuon spectra to medium-induced effects, we consider two extreme cases: 
(i) the "minimum" effect with collisional energy loss only and (ii) the "maximum" effect with collisional 
and radiative energy loss in the incoherent limit of independent emissions without taking into account 
the LPM coherent suppression of radiation (i.e. dE/dx cx E and is independent of path length, L) ]l3| . 
In the latter scenario we use the Bethe-Heitler cross section obtained in relativistic kinematics and derive 
the medium-induced radiative energy loss per unit length of quark of mass M q [fL9| as the integral over 
the gluon radiation spectrum 
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l-MJE 

_ p f , 4g s C 3 (y)(4-4; / + 2y 2 ) 
dx ^ P J aV 9ny[My + ml(l-y)Y {) 

9na 2 C ab f 2 2] 2(g^(l-y)) 1/4 
C 3 (y) = - [1 + (1 - y) - y \ In , 

where m g ~ 3T is the effective mass of the emitted gluon at temperature T; ?/ is the fraction of the initial 
quark energy carried by the emitted gluon; p oc T 3 is the density of the medium; C ab = 9/4, 1 and 4/9 
for and qq scatterings respectively; a s is the strong coupling constant for Nf active quark flavours; 

the Debye screening mass p,o regularizes the integrated parton rescattering cross section. 

The dominant contribution to the differential cross section da/dt for scattering of a quark with energy 



E and momentum p = y E 2 — M 2 off the "thermal" partons with energy (or effective mass) mo(r 
3T(r) <C E at temperature T can be written in the target frame as [Q 

da ab _ „ 27ra 2 (t) £ 2 _ ^ _ 12tt 

f+ / /, 

1 QCD 



dt ~ Cab t> p 2 ' a s{t) — ^33 _ 27V/) In (£/Aq CD ) ' (4) 



and QCD scale parameter Aqcd on the order of the critical temperature, Aqcd — T c . The integrated 
parton scattering cross section has the form: 

*a b {T) = j dt^ (5) 

where t max = [s - (M q + m ) 2 ][s - (M q - m ) 2 ]/s and s = 2m E + to 2 , + Af 2 . 

In the z-th rescattering off a comoving medium constituent (i.e. with the same longitudinal rapidity y) 
with squared momentum transfer U and effective mass m$i , the quark loses total transverse energy Aeri 
and change rapidity on Ay: 
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kt cos < 



Et — Aer 

where the transverse momentum kick per scattering is 



— TT--£-) +Afc 2 sin 2 ^ + M 2 , (6) 
p T 2m 0l 2p T J 

(7) 




(8) 

and the angle <p between the direction of vector k t i and axis z being distributed uniformly. The medium- 
induced radiative energy loss is calculated with Eq. (^|) without modification of the longitudinal rapidity. 

In our calculations, we use the Bjorken scaling solution Jl7|] for the space-time evolution of the energy 
density, temperature and density of the plasma: 

£ (r)r 4 / 3 = e r 4/3 , T{t)t^ = T r^\ p(r)r = p r Q . (9) 

To be specific, we use the initial conditions for a gluon-dominated plasma expected for central Pb+Pb 
collisions at LHC Q: t ~ 0.1 fm/c, T ~ 1 GeV, N f w 0, p g « 1.95T 3 . It is interesting that the 
initial energy density, Eq, in the dense zone depends on b very slightly, feo/eo 10%, up to b ~ Ra 
and decreases rapidly for b > Ra [§0|- ( - >n * ne °th er hand, the proper time of a jet to escape the dense 
zone averaged over all possible jet production vertices, (tl), is found to decrease almost linearly with 
increasing impact parameter. This means that for impact parameters b < Ra, where 60% of the heavy 
quark pairs are produced, the difference in rescattering intensity and the corresponding energy loss is 
determined mainly by the different path lengths rather than the initial energy density. 



Finally, the nuclear shadowing corrections according with EKS model |21| have been also taken into 
account in our calculations. 

The simulation of quark rescattering is halted if one of the following three conditions is fulfilled: 
1) A quark escapes from the dense zone, i.e. its path length becomes greater than the effective transverse 
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spread of the matter from the production vertex to the escape point. The details of the geometrical 
calculations of these quantities at a given impact parameter can be found in Ref . [^0| . 

2) The plasma cools down to T c — 200 MeV. We thus neglect possible additional small contributions to 
the total energy loss due to re-interactions in the hadron gas. 

3) A quark loses so much energy that its transverse momentum px drops below the average transverse 
momentum of the "thermal" constituents of the medium. In this case, such a quark is considered to be 
"thcrmalized" and its momentum in the rest frame of the fluid is generated from the random "thermal" 
distribution, dN/d 3 p oc exp (-E/T), boosted to the center-of-mass of the nucleus-nucleus collision fll| , |l2|| . 

Figure 3 shows the invariant mass spectra from BB decays for various nuclear effect scenarios: 

(1) no loss and shadowing; (2) no loss, with shadowing; (3) with collisional loss and shadowing; (4) with 
collisional and radiative loss and shadowing. The absolute normalization is the same as in previous section. 
The results for all quark production mechanisms are qualitatively similar: the shadowing corrections are 
relatively small (at the level of « 15%); the collisional loss reduces dimuon rate by factor ~ 1.3 — 1.6; the 
additional radiative loss reduce the rate up to factor ~ 3-4. The relative contribution of radiative loss 
grows with increasing M and pr due to the stronger energy dependence of the loss. The dimuon spectra 
from direct quarks are some more sensitive to the energy loss of 6-quarks as compared with dimuons from 
"showering" quarks due to different kinematics (e.g. strong azimuthal back-to-back correlation). 
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FIG. 3. Invariant mass distribution of pairs 
from BB decays for various scenarios: without loss 
and shadowing (dash-dotted), without loss and with 
shadowing (solid), with collisional loss and shadowing 
(dashed), with radiative and collisional loss and shad- 
owing (dotted). From top to bottom: "pair creation", 
"flavour excitation" and "gluon splitting". 




10 12 14 16 , 18 20 22 24 



p„ GeV/c 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



FIG. 4. Summarized over all bb production 
mechanisms (a) transverse momentum and (b) 
pseudo-rapidity distributions of B — > J/ip(—> /i + /i~) 
decays for various scenarios: without loss and shad- 
owing (dash-dotted), without loss and with shadowing 
(solid), with collisional loss and shadowing (dashed), 
with radiative and collisional loss and shadowing (dot- 
ted). 



The summarized over all bb production mechanisms transverse momentum and pseudo-rapidity dis- 
tributions of B — > J/tp(— > /x + /i~) decays is presented in figure 4. Including nuclear shadowing reduces 
secondary charmonium yield by ~ 25%, while the final state rescattering and energy loss by 6-quarks can 
further reduce the J/ip rates by a factor of ~ 1.3-2 in the CMS kinematical acceptance. 

We see that the influence of nuclear effects on secondary J/ip's and high-mass dimuons are quite 
different: the decrease of the B — > J/ip rate by nuclear shadowing is comparable to the effect of medium- 
induced final state interactions. The increased sensitivity to nuclear shadowing is due to the different x 
and Q 2 regions probed. The different influence of energy loss on secondary charmonia and dimuons from 
BB decays is because J/i/Vs come from the decay of a single b quark instead of a bb pair and there is a 
non-negligible probability that the energy lost by one quark is small. Thus since secondary charmonia 
reflect energy loss of only one 6-quark, the corresponding suppression is less than for semileptonic BB 
decays and is independent on quark production mechanism. Thus we believe that a comparison between 
high-mass dimuon and secondary J ftp production could clarify the nature of energy loss. 
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IV. CONCLUSIONS 



We have considered the different mechanisms of heavy quark production in ultrarelativistic nuclear 
collisions and found that the " showering" bb pairs at LHC give the same order of contribution to high-mass 
dimuon spectra (Mu+u- > 20 GeV/c 2 ) as direct quark pairs. The bb pair production from parton showers 
is dominant mechanism for B — > J/t/j production. The medium-induced parton rescattering and collisional 
energy loss can reduce high-mass dimuon rate by factor ~ 1.5, and the additional radiative energy loss can 
reduce the rate up to factor ~ 3-4 (depending on quark production mechanism) . The relative contribution 
of radiative loss grows with increasing M and pr due to the stronger energy dependence of the loss. Due 
to different kinematics the high-mass dimuon spectra from direct quarks are more sensitive to the energy 
loss of 6-quarks than dimuons from "showering" quarks. The nuclear shadowing corrections is w 15% for 
high-mass dimuons. Since secondary charmonium production reflects the energy loss of only one 6-quark, 
the corresponding suppression by a factor of 1.3 — 2 is lesser than for BB decays and is independent on 
quark production mechanism, but shadowing corrections are some larger. 

We conclude that the dimuon spectra will be sensitive to final state rescattering and energy loss of 
bottom quarks in dense matter. However, there are still theoretical uncertainties in the initial production 
of heavy flavours in nucleon-nucleon collisions at LHC energies. Thus measurements in pp or dd collisions 
at the same or similar energies per nucleon as in the heavy ion runs are required. 
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